Boundary effects in the density-matrix renormalization group calculation 
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We investigate the boundary effect of the density matrix renormalization group calculation (DMRG), which 
is an artifactual induction of symmetry-breaking pseudo-long-range order and takes place when the long range 
quantum fluctuation cannot be properly included in the variational wave function due to numerical limitation. 
The open boundary condition often used in DMRG suffers from the boundary effect the most severely, which 
is directly reflected in the distinct spatial modulations of the local physical quantity. By contrast, the other 
boundary conditions such as the periodic one or the sin 2 -deformed interaction [A. Gendiar, R. Krcmar, and 
T. Nishino, Prog. Theor. Phys. 122, 953 (2009)] keep spatial homogeneity, and are relatively free from the 
boundary effect. By comparing the numerical results of those various boundary conditions, we show that the 
open boundary condition sometimes gives unreliable results even after the finite-size scaling. We conclude that 
the examination of the boundary condition dependence is required besides the usual treatment based on the 
system size or accuracy dependence in cases where the long-range quantum fluctuation is important. 

PACS numbers: 71.10.Hf, 71.27,+a, 71.10.-w 
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I. INTRODUCTION 

The density matrix renormalization group (DMRG) pro- 
posed by Whitei is by now one of the most powerful nu- 
merical tools to determine the ground states and low-energy 
excitations 2 ^ of the correlated quantum many-body sys- 
tems in low dimensions. Accumulated studies revealed that 
the DMRG has advantages when the system is not critical 
(gapped), where one could effectively reduce the number of 
bases states without sacrificing the numerical accuracy. The 
open boundary condition (OBC) usually used in the DMRG 
analysis also works to reduce the number of basis, which 
was pointed out by White in the first staged In fact, com- 
pared to the periodic boundary condition (PBC), OBC has a 
smaller number of joint points between the blocks (subsys- 
tems) than PBC (one and two points, respectively), and thus 
has smaller entanglement entropy which follows a concept of 
area law»i the smaller the size of the boundary between the 
two blocks, the smaller the entanglement entropy one finds 
between them. However, one shall have to be careful in this 
context because the existence of the open edge points does not 
only indicate the reduced entanglement, but means the break- 
ing of the translational symmetry which leads to two artifacts 
in critical (gapless) systems. One is to favor a particular frac- 
tion of ground-state manifolds which are equivalently degen- 
erate in the bulk limit. The other is to mix partially the states 
which have different symmetric properties from the ground- 
state manifold. The former effect does not modify the ground- 
state nature (only chooses a pure state) and is equivalent to 
introducing the infinitesimal symmetry-breaking field in large 
systems (see Sec. II A for details). However, the latter may 
sometimes allow some low-energy excited state to overwhelm 
the original ground state. For example, when the correlation 
length is longer than the typical system size, the oscillation 
induced by the open boundary may open an apparent gap, par- 
ticularly when the number of states kept is not large enough. 
This is what we call the "boundary effect." 

The portion of boundary is order 1 while the system is or- 



der L, and by the finite-size scaling, the boundary effect is 
expected to decay faster than 1/L on an average. Therefore, 
by the proper finite-size scaling, one could in principle get rid 
of the boundary effect, and the correlation with the largest 
length scale is considered to be the one characterizing the 
ground state. In fact, one can take advantage of the OBC in 
order to examine the nature of orders or correlations: Due to 
the boundary edges, the wave function loses the translational 
symmetry and the two-point correlation of the largest length 
scale explicitly appears as a spatial modulation of local phys- 
ical quantities. This treatment was actually applied to DMRG 
by the present authors. In the one-dimensional Kondo lattice 
model (KLM) the two open edges are regarded as impurities 
which induce the Friedel oscillations, and by measuring the 
structural factors of the charge or spin amplitudes, the authors 
detected the wave numbers, 2Uf or of the Tomonaga- 
Luttinger liquid state. 5 Also in two-dimension, the competi- 
tion of several different types of long range orders with differ- 
ent spatial periodicity is detected by the spatial modulation of 
particle density in the fermionic model on an anisotropic tri- 
angular lattice under the chemical potential on system edges^ 
Such analysis is, however, applicable only when the basic na- 
ture of the ground state is already understood; namely, it clar- 
ifies the details of the already known (stable) ground state but 
not the ground state itself. Moreover, the analysis requires a 
set of results with high enough accuracy as well as a parame- 
ter region where the scaling law is safely adopted. Once either 
of the above two is not fulfilled, one may no longer obtain a 
reliable result, which is often the case in DMRG, since there 
is an upper bound of the number of states kept, in practice. 

The present paper discusses how to classify systematically 
the state which does or does not break any symmetry of 
the original Hamiltonian within a usual numerical accuracy 
of DMRG. We examine how the critical behavior of quan- 
tum many-body states at finite system size is influenced by 
the boundary condition by choosing two characteristic exam- 
ples: a well-known symmetry-broken long-range order in a 
spin chain, and strongly correlated electronic states with ex- 
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tremely long correlation length. We consider the modified 
or deformed OBC, as well as PBC, and demonstrate that the 
expectation value of the local operator depends severely on 
boundary conditions so that they cannot be a reliable mea- 
sure. We finally see that in order to conclude the presence of 
long-range order, it is necessary to confirm that the two-point 
correlations remain finite toward the bulk limit irrespective of 
the boundary conditions. 

II. BOUNDARY CONDITIONS 
A. Symmetry-breaking long-range orders 

We consider a class of long-range orders due to transla- 
tional symmetry breaking; i.e., the order operators do not 
commute with the quantum Hamiltonian of crystals which 
usually keeps the translational symmetry. It is known that 
even when such spontaneous symmetry breaking occurs in 
the bulk limit, quantum fluctuation "obscures" the breaking of 
symmetry at finite system size and one finds a unique ground 
state with perfect symmetry^ 

Let us briefly follow the contexts of Ref. by Koma and 
Tasaki, which is based on two different kinds of measures, 

m s = lim lim m s (B,L), 

rn s (B,L) = ^(Y,Oi{B,L)), (1) 

i 

a = lim 

L— foo 

°{L) = ^((E a '( B = ' L )) 2 )' ( 2 ) 

where (Oi(B, L)) is the expectation value (trace) of local op- 
erator Oi(B, L) on the ith site under the symmetry-breaking 
field B and for system size L. Hereafter, we use mere O, as 
the one with B = and finite L except otherwise noted. In the 
Heisenberg antiferromagnet, m s is the staggered magnetiza- 
tion obtained by applying the infinitesimally small symmetry- 
breaking field, B — >• 0, whereas the latter is the conventional 
long-range order parameter which consists of two-point cor- 
relation functions, (0,Oj), between sites i and j. Koma and 
Tasaki proved that 

m s > ra, (3) 

holds for constant r at all temperature for various types of 
models with long-range order [r = 1 holds in general, but a 
higher symmetry gives stronger bounds, e.g., r = \/3 for the 
SU(2) Heisenberg antiferromagnet^]. Here, a > means 
that the two-point correlation function does not decay to zero 
in the bulk limit, (OiOj) ^ for \i — j\ — > oo. Therefore, 
Eq.([3j indicates that a > guarantees the existence of sym- 
metry breaking represented by m s > 0. 

We now interpret Eq.© to the numerical analysis on a fi- 
nite size system at zero temperature. Note that the following 
Eqs. (|4]i-(|7]i and conditions (I) — (III) are conjectures, which are 



derived logically in the following part of this section. Since 
the finiteness of the system is characterized by the presence of 
boundaries, we first discuss the role of boundary conditions. 
When the system does not have any open edges, e.g., peri- 
odic (PBC) or antiperiodic (APBC) boundaries, the transla- 
tional symmetry in the Hamiltonian "obscures" the breaking 
of translational symmetry at finite L, and we find (Oi(L < 
oo))pbc = 0. Thus we need to include the infinitesimal 
symmetry-breaking field B, or analyze the two-point corre- 
lation function o-(L) pbc , which coincides with a in the limit 
of large L: lim a(L) phc = a. 

L— >oo 

In the case of open boundary condition, the coupling be- 
tween the 1st and the iVth sites in PBC is missing and the 
translational symmetry is broken in the Hamiltonian. This 
works to discriminate the two edge bonds energetically from 
the rest of the bonds in the system, which approximately cor- 
responds to placing the effective external field on the edge 
bonds Such effective local field induces oscillation of local 
correlations, which is absent in PBC. 

We basically confine ourselves to the case where the pe- 
riod of the oscillation is compatible with the interval between 
edge bonds; e.g., if the induced oscillation has twofold peri- 
odicity, the number of bonds, L — 1, must be odd. Then the 
open boundary condition corresponds to a locally introduced 
symmetry-breaking "field," B ob , which is coupled to the lo- 
cal operator Oi = o.; + io.; + 2 — o.i&i+i at the edge sites, i = 1 
and N, where o; is the operator acting on the ith site. The 
effect of B b decays as it propagates toward the system center 
from both edges, which is reflected in the gradual decrease of 
the oscillation amplitude of (Oj) with i N/2. Then, we 
expect m s (L) bc > m s . The two-point correlation functions 
(OiOj) hc are those under the effective symmetry breaking 
"field," so that we shall find u(£) bc > cr. 

Basically the portion of boundary is order 1 while the sys- 
tem is order L, and taking the limit of L — > oo corresponds 
qualitatively to having B \> —> as in Eq.([T]i. Thus both of 
the above two boundary conditions (PBC and OBC) lead to 
the same conclusion, 

m s = lim lim m s (B,L) pbc (4) 

B— >0 L— >oo 

= lim m s (B = 0, L) ohc , (5) 

L—>OG 

a = lim a(L) phc (6) 

L— yoo 

= lim o-(L) ob c, (7) 

L— f oo 

with m s > a. We therefore conclude that when a > 0, we 
always find m s > 0. 

However, the numerical approximation, e.g., restricting the 
number of basis to m in DMRG, sacrifices the inclusion of 
long-range quantum fluctuation to some extent. The OBC 
does not only work as an "effective symmetry breaking field," 
but its locality induces a nonuniformity of the quantum fluc- 
tuation and mixes the states which have different symmetry 
from the original pure state. While such additional effects 
may be weakened as L — > oo, incomplete long-range quan- 
tum fluctuations sometimes cause harmful effect on the origi- 
nal ground-state symmetry as we will see shortly. Therefore, 
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the practical procedure to systematically get rid of the artifact 
of the boundary condition is required. 

The above-mentioned nonuniformity under OBC is re- 
flected in {Oil and {O n O n +i), and their deviation from the 
one in the uniform system is enhanced near the open edges. 
Since we need to analyze the systematic behavior of these 
quantities in large systems, we examine them at each site i and 
n instead of taking their mean value over the whole system as 
in Eqs. ([T} and (|2}. Then, those of site i and n away from 
the edges are adopted which should fulfill lirnj^o^Oj) > 
rlim L ^ 00 (0„O n+i ) as in Eq. ©. 

With this in mind, we propose the necessary conditions the 
numerical results must fulfill in order to safely conclude the 
breaking of translational symmetry in the bulk limit: 

(I) the two-point correlation function {O n O n +i) remains fi- 

nite for i — > L/2, after the L and m scaling, 

(II) the conclusion obtained by (I) does not depend on the 

boundary conditions, including those with and without 
symmetry breaking, 

(III) as {Oi) of finite systems becomes uniform under the 

variation of the boundary conditions, the corresponding 
two-point correlation function {O n O n +i) approaches 
the value of the pure state in the bulk limit. 

The rest of the paper is devoted to the numerical "proof" of 
this conjecture or proposal in the representative models in one 
dimension. 

The boundary conditions we consider are classified into 
those which keep the translational symmetry and those which 
do not. As a typical boundary condition of the former class, 
we deal with the system with "deformed interactions" recently 
proposed by Gendiar, Krcmar, and Nishinor^ as well as the 
PBC and APBC. The representative condition of the latter 
class is the OBC. Besides, we consider the cases referred to 
as modified open boundary, which is the open boundary with 
potentials on edge sites or edge bonds first adopted by White 
and Huse.— As we discussed in this section, the open bound- 
ary works as an effective field on edges, which brings the spa- 
cial nonuniformity of the quantum states. This artificial effect 
can be suppressed by adjusting the bond strength or potentials 
on edges by hand, and thus one can tune the degree of inho- 
mogeneity by the "modification of the open boundary." 

We finally comment on the cases where the period of oscil- 
lation of local quantity {Oi) does not match the length of the 
system L. For example, in the case of dimer order, the local 
operator Oi is defined as Oi = Sj+i • Si+2 — Si • S»+i, and the 
twofold periodic oscillation appears in the nearest-neighbor 
spin-spin correlation. In the odd-L system with PBC, a kink, 
namely a twisting oscillation emerges in the system. Since the 
PBC keeps translational symmetry, the wave function should 
be the superposition of L different wave functions, each with 
a kink on the ith (i = 1 ~ L) bond. Then, {Oi) is uniform 
but is suppressed from the one without a kink by the order of 
1/L, which recovers the bulk value when L^co. 

When the system of odd L has an open boundary, this kink 
is confined to the center site as (O^/a) ~ 0. This is because 
the dimer bond is pinned the strongest at both edge bonds, and 



the oscillations of (0,) starting from these strong edge bonds 
interfere at the center and form a kink. — The effect of this kink 
on {Ol/2) cannot be excluded by the m or L scaling, since it 
is locked at the center of the system. The artifact of kink on the 
value of (O n O n +i) is also nontrivial, and less easy to get rid 
of compared to the usual boundary effect of OBC with even L. 
We show in Sec. IV that {O n O n +i) in the dimer state with a 
kink in OBC finally approaches the value without a kink with 
phase shift tt, which means that (II) basically holds regardless 
of the presence or absence of a kink. In this sense we can 
judge the presence of the dimer order even in the system of 
odd L with a kinkJi 

In the above semiclassical picture, the kink is localized on a 
single bond. However, in the quantum systems we deal with, 
the kink spreads over a certain length scale due to quantum 
fluctuation, and this length scale depends on the stiffness of 
the dimer order. Still, if we take L enough larger than this 
lengthscale, the same discussion holds. 



B. Modified open boundary 

The analysis with "modified open boundary" in DMRG was 
first adopted in the 5 = 1 Haldane chain with 5=1/2 spins 
on the edge sites by White and Huse^ To evaluate the Hal- 
dane gap of infinite systems, they adjusted the coupling, J e nd> 
between the 5=1/2 edge spin and the neighboring 5=1 
spin, and minimized the nonuniformity of the ground state 
away from the edge. 

In the present paper, we start from the usual OBC and mod- 
ify the amplitude of bond interactions on both ends or place 
the potential on both edge sites to analyze the boundary con- 
dition dependence of the ground state. As for the spin system, 
the bonds, J c dgc, on both left and right edges are varied. If 
we place the strong antiferromagnetic bond on one edge, the 
singlet correlation on that edge is enhanced, and as a result, 
the correlation on the neighboring bond decreases, which en- 
hances or reduces the inhomogeneity of the ground state. In a 
similar manner, the electronic states can be tuned by modify- 
ing the hopping energy of edge bonds, t c ^ gc , or by placing the 
chemical potential on edge sites, /x e dge- One can minimize the 
nonuniformity of the electron density by adjusting the value 
of /i dgc- Furthermore if one uses different values of /i c d gc 
between the two edges, the conduction electrons shift to left 
or right, and the position of the center of mass of electrons can 
be controlled. 

Besides the system with even number of L, we deal with 
odd number of L to shift the center of mass of electrons 
smoothly along the ID chain in OBC and at the same time 
keeping electron density constant away from the boundary. 
Although the average electron filling factor p = N e /L (N e 
denotes the electron number) slightly deviates from the one 
we have in the even-L case, we can pin the excess electrons 
or holes to the boundary by tuning the value of ^ e dge- 
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FIG. 1: (a) Site dependence of the nearest-neighbor correlation, 
(SiSi +1 ), in the 5 = 1/2 Heisenberg model with OBC and sin 2 
deformation at L = 32 and m — 200. (b) The functional form of 
the sin 2 deformation in Eq.([8]l, and the two types of wave function 
of the dimer state, ipi and tpz, which have the same energy. 



C. Deformed interactions 

Recently, Gendiar, Krcmar, and Nishino proposed an un- 
precedented analysis to get rid of the boundary effect by de- 
forming the interaction strength of the system following the 
sin 2 function which decreases from the center toward both 
ends of the system . 11 ' 16- — 

In this paper we use the following deformed Hamiltonian: 



L-l 



H = ^ sin 2 (— jg(n,n+l) 

n=l 

^ . 2 /7r(n + l/2)\ t . 

n=l ^ ' 

A . 2 /7r(n-l/2)\ , , 



L 



(8) 



where g(n, n+1) is the nearest-neighbor interaction, h(n, n+ 
2) is the next-nearest-neighbor interaction, and u(n) includes 
the on-site interaction and potential. The uniform chemical 
potential, /i, is also included as u(n) which is deformed by 
the prefactor and becomes nonuniform^ Figure Q] shows the 
site dependence of the nearest-neighbor spin-spin correlation 
of the 5* = 1/2 Heisenberg model, g(n, n + 1) = J S„-S„+i, 
under the usual OBC and under the deformation of interaction. 
One finds that the usual OBC induces a twofold oscillation 
which decays slowly toward the system center. The oscillation 
is fully suppressed by the deformation and the local quantity 
becomes site independent. 
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FIG. 2: Comparison of local correlation ,{Sj Z Sf +1 ), in the J1-J2 
model calculated under usual OBC, modified OBC, deformed OBC, 
and PBC, with L = 32 and m = 200. The convergence of the results 
regarding the m dependence is confirmed. For modified OBC we use 
Jcdgc = 0.475J on both ends. The results of OBC with L = 96 is 
added to compare with the L = 32 one. 



Here, we give an interpretation on these "site-independent" 
results. Consider the system with even L = 2N. Then, one 
finds the following relation: 



N 

E 

1=1 



J sin' 



n(2l) 



N 



1=1 



tt(2Z - 1) 



(9) 



which means that the sum of the coupling constant is the 
same between even bonds (21) and odd bonds (21 — 1) [al- 
though the total number of even bonds with finite coupling 
is less than that of odd bonds because sin 2 (27V7r/L) = 0, 
which corresponds to the missing bond in OBC]. Suppose 
that we have translational symmetry broken dimer states with 
twofold periodicity, which are represented by two different 
types of wave functions, ipi and ipi as shown in the lower 
panel of Fig. [TJ 991 consists of dimers on even bonds and 
<f 2 consists of those on odd bonds. If either of ipi and (p2 has 
lower energy than the other, that state is selected as the ground 
state. However, due to the relation in Eq.fO, tpi and ip2 have 
the same energy, and thus the translationally symmetric wave 
function (tpi ± tf2)/V2 can be constructed, which has the 
same ground-state energy. Hikihara and Nishino recently cal- 
culated the electronic system with this sin 2 deformation and 
found that the overlap with the PBC wave function is almost 
1 M This result was in fact supported analytically for the non- 
interacting case. 21 In the results of free fermionic modelr^ 
the Heisenberg chain [Fig. [He)], and 3\-Ji models, which we 
will see shortly, the translational symmetry is indeed recov- 
ered by the deformation. Therefore, we consider that the sin 2 
deformation is a condition to have the translationally symmet- 
ric ground state at least for simple single-band short-range in- 
teracting systems. 
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FIG. 3: Comparison of two-point correlation functions, 
\(0„O n +i)\, in the J1-J2 model with L = 196 and m = 300. 
Those of usual OBC, modified OBC, and deformed OBC start from 
n = 40, and the PBC one from n = 1. The convergence of the 
results regarding the m dependence is confirmed. Inset shows the 
results of modified OBC with n — 60 and 40. 



m. J1-J2 MODEL 
A. Dimer order 

Our first step is to examine how the order operators dis- 
cussed in Sec. Ill Al behave under the variation of boundary 
conditions when the symmetry-breaking long-range order is 
present in the bulk limit. For this purpose, we choose a J\- 3% 
model, whose Hamiltonian reads 

U = s £ j (J x s j - s i+i + J 2 Sj-S i+2 ) , (10) 

3 



B. Boundary condition dependence 

We now compare the results under three different types of 
boundary conditions at J2/J1 = 0.4 where the dimer long- 
range order is present in the bulk limit. Figure [2] shows the 
site dependence of the local operator, (S?Sf +1 ) = (Di)/3. 
Under usual OBC, (Di) shows twofold periodic oscillation 
with a large amplitude. The amplitude of oscillation is signif- 
icantly suppressed, when the modified OBC is adopted with 
^edgo = 0.475J on both edge bonds. Further, by the sin 2 
deformation, the oscillation is almost completely suppressed 
and the spatially uniform (SfSf +1 ) is obtained, which means 
that the translational symmetry is recovered. In fact, this de- 
formation gives almost the identical results with that of the 
PBC given in the same figure. In all cases, the magnitude 
of I (5? S? +1 ) I has only small i dependence. We also plot 
the results under OBC with L = 96 to show that the L de- 
pendence of the oscillation amplitude is significantly small, 
indicating that the distinct twofold oscillation remains after 
the finite-size scaling in OBC. Therefore, there is no doubt 
that the extrapolated value of (SfS? +1 ) in the bulk limit, 
liniL^oo(-DL/2). depends severely on the boundary condi- 
tions. 

Next, we show two point correlation functions, 
\(O n O n +i)\, in Fig. [3] under the same choices of boundary 
conditions used in Fig. [2] One finds that all cases almost 
asymptotically approach a constant value after i > 30, in 
sharp contrast to the severely boundary-dependent (Di). The 
slight difference for the modified boundary condition is due 
to a boundary effect of finite systems, and it almost vanishes 
for large n > 60 as shown in the inset of Fig. [3] From this 
result, one can safely confirm that the extrapolated value of 
the two-point correlation function, iimi_ ) . 00 \(O n O n +i)\, is a 
reliable measure of symmetry-breaking long-range order for 
any choices of boundary conditions. 



where J\ > 0. The exact ground state at a = J2/J1 = 0.5, 
which is called the Majumdar-Ghosh (MG) point, is exactly 
represented by the product of local dimer states, and thus 
the dimer-dimer correlation is expected to show no decay for 
length scale large enough to neglect the boundary effect in- 
duced by OBC. The existence of dimer long-range order is 
established off the MG point in the region of a > 0.24:^- 
The nearest-neighbor spin-spin correlation operator is given 

as 

A = S 4 -S m , (11) 

where S, is the spin operator of site i. We use the dimer oper- 
ator 

Oi = Di + \ — Di = Sj + i • S; + 2 — Sj- Sj + i (12) 

to define the order parameter a in Eq. ©. When the trans- 
lational symmetry breaking dimer order is present, the two- 
point correlation function, (OiOj), remains finite in the bulk 
limit. 



IV. KONDO LATTICE MODEL 

A. Preliminary information 

The next example is devoted to the Kondo lattice model, 
which is one of the basic models for the heavy-fermionic 
systems 2 - and is studied also as a prototype system to clar- 
ify the effects of coupling of localized spins and conduction 
electrons j2i The Hamiltonian in one dimension is given as 

n = t E ( c *rc;+i* + H - c -) • 

3 j 

where Cj a denotes the annihilation operator of the conduction 
electron at the jth site with spin a =| , I, and Sj denotes the 
localized spin operator with 5 = 1/2. Sj = ^T^^icy-Cja' is 
the spin operator of the conduction electron with Pauli matrix 
T a tr i. The conduction electrons hop with energy t and in- 
teract with localized spins through the Kondo exchange cou- 
pling J. Its ground state in one dimension is considered as 
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a Tomonaga-Luttinger liquid (TLL) in the weak-coupling re- 
gion away from half filling. In the strong-coupling region, 
each conduction electron forms a singlet with localized spins 
and behaves as a single hole, and the ferromagnetic metallic 
ground state is realized^ 

A few years ago, however, there was a proposal by Xavier 
et.al. on the possible dimer phase in the region of J/t < 1.6 
at quarter filling,™ which was originally considered a param- 
agnetic TLL region. Based on the fact that (Ol/2) calculated 
under OBC with DMRG seems to remain finite after the size 
scaling, they claimed that the translational symmetry is broken 
in the bulk limit. We carried out similar DMRG study on the 
same model and reproduced their results, but with different 
conclusions based on the analyses of the boundary effecti^l2& 
In fact, the model at < J/t < 1 is in a state which is ex- 
tremely difficult to analyze even by DMRG established in one 
dimension^ 

In the following we show in detail how the two-point cor- 
relations as well as the local quantities behave under the vari- 
ation of boundary conditions, which has distinct differences 
from those of the dimer phase in the J\-Ji model. 



B. Modified open boundary 

In the KLM the open boundary condition can be modified 
by the chemical potential /i d g0 = m(1 — $) and fj,6 on the left 
and the right edges, respectively. The value of fi/t is tuned 
to —1.6 in order to keep the density of charges at quarter fill- 
ing within the accuracy of 10 -3 away from the edge sites, 
and 5 is introduced to shift the center of mass of electrons. 
Figures Ua) and Ub) show the site dependences of the bond 
kinetic energy (cjcj_|_i + h.c.) and the nearest-neighbor corre- 
lation of localized spins (S^S^ +l ) for several choices of 8. 

Under the usual OBC with even L, both (cfcj+i + h.c.) 
and (5f5f +1 ) show large oscillation. The same oscillation 
is found when the modified boundary is adapted to the odd- 
L system with 6 = 0, namely when /.t is placed on only 
one of the edge sites. From this result, one can under- 
stand that fi works to pin the extra charges off quarter filling, 
N e — (L — l)/2, to the edge site, and the rest of the system 
remains quarter filled. Then, as 6 is increased from toward 
0.5, the amplitude of oscillation is gradually suppressed. 

Figure Uc) shows the schematic illustration of the effect of 
/i and <5 on edge sites. The usual OBC with an even number of 
lattice sites pins the electron density to be the bond-centered 
spatial modulation along the one-dimensional chain. When 
the boundary is modified as < 5 < 0.5, the electrons on the 
left and right edges have different hopping amplitude to their 
neighbors, which works to displace the spatial modulation of 
hoppings and to suppress the oscillation of hopping amplitude 
as we find in Fig. EJa). Since the hopping of electrons in- 
duces correlations between the local spins through the Kondo 
exchange interaction, 5 also modifies the local spin correla- 
tions as shown in Fig.HJb). This result arises the question of 
whether the oscillation of (S*S* +1 ) found in the OBC is an 
intrinsic property of the KLM. To study the boundary con- 
dition dependence in detail, we next analyze the two-point 
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FIG. 4: Spatial modulation of two local quantities, (a) bond kinetic 
energy, (c}ci+i + h.c), and (b) nearest-neighbor correlation of lo- 
calized spins, {SiSf + i}, in the Kondo lattice model calculated with 
usual OBC with L = 24 and modified OBC with L = 25 with 
m = 600. The convergence of the results regarding the m depen- 
dence is confirmed. For modified OBC we use /i e d gc = m(1 — 5) 
and fiS on left and right edge site, respectively, with fj,/t = —1.6. (c) 
Schematic illustration of the electronic states under modified bound- 
ary condition by the chemical potential, /i c d g c • The size of the shaded 
circles represents the density of conduction electrons. 



correlation function of the local spins. Figure [5] shows the 
dimer-dimer correlation of the local spins \(O n O n +i)\. For 
usual OBC, \(O n O n +i) \ saturates toward the finite value al- 
ready at i ~ 20, which apparently suggests the existence of a 
dimer state. However, when the boundary is modified (5 > 0), 
the functional form changes significantly and starts to decay 
monotonously. 

To examine whether the odd L of the one-dimensional 
chain artificially suppresses the dimer correlation or not, we 
also calculate | (O n O n +i) \ of the J1-J2 model under both even 
and odd L starting from n = 5 as in KLM, which is shown 
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FIG. 5: Two-point correlation functions of the dimer operator, 
\{O n On+i)\ (for fixed n — 5) in the Kondo lattice model with 
m — 600 (the m dependence is found to be negligible for the present 
L's). The usual OBC with L = 32 and modified OBC with L = 33 
for several choices of 5 (/Wge is the same as Fig|4j are compared. 
The inset shows the two correlation functions, \{O n O n +i)\, of the 
J\-Ji model with usual OBC at L — 32 and 33 with n = 5, to be 
compared with those of the Kondo lattice model in the main panel. 



in the inset of Fig|5] In the case of the J\-Ji model with the 
true long-range dimer order in the bulk limit, | (O n O n +i) | of 
the odd L first decays and takes the minimum (or a dip) at 
around i ~ 12 which is near the system center, then shows an 
upturn, and finally approaches the value of the even L, which 
is almost an i-independent constant. The local decay (dip) of 
the dimer correlation is due to the kink structure formed by 
the phase shift 7r between the two dimer orders starting from 
both ends: the misfit of odd L and the twofold periodicity of 
the dimer order. Since the kink is located around the center of 
the system, the dimer correlation gradually decreases toward 
the system center. But for long distances beyond the center of 
the system, this kink is always located between the two dimer 
operators Oi and Oj, and the dimer correlation increases up 
to its original value with an additional phase shift by tt. 

Such nonmonotonic correlation function is also expected 
for the KLM of odd L, if the long range dimer order is present. 
However, the results in Fig|5]for finite 6 show that the dimer 
correlation monotonically decreases for long distance even for 
the case of symmetric edge potentials at 6 — 0.5. This result 
suggests that the dimer order is not an essential property of the 
KLM at quarter filling at J/t = 1. The large <5 dependence of 
electronic states [see Fig. Ufa)] and concurrently of the spa- 
tial structure of localized spins [Fig. SJb)] are the sign that 
the electronic state is severely boundary condition dependent, 
and one should analyze the finite-size scaling and its boundary 
condition dependence in order to determine the bulk property 
of the system safely. 



C. sin deformation and periodic boundary 

The final set of results is devoted to examining the effect 
of finite size L and finite number of basis states m in DMRG 



FIG. 6: Site dependence of the electron density in the quarter-filled 
Kondo lattice model at J/t = 1, under the usual OBC and sin 2 
deformed OBC with L = 32 and m = 600. 
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FIG. 7: Site dependence of (a) bond kinetic energy, (ejci+i + h.c), 
and (b) nearest-neighbor correlation of localized spins, (5j S^^i), in 
the quarter-filled Kondo lattice model at J/t — 1, under the usual 
OBC and sin 2 deformed OBC with L = 32 and m = 600. 



under the condition where the effect of the missing bond in 
OBC is as weakened as possible. Based on the discussion 
in Sec. Ill CI we perform the calculation under the sin 2 defor- 
mation, which reproduces the translationally invariant wave 
function for the J1-J2 model. 

In the KLM, we add uniform chemical potential as n- 
independent u(n) which is needed to keep the electron den- 
sity at quarter filling, and deform the Hamiltonian according 
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FIG. 8: Two-point correlation functions of the dimer operator, 
\(0„O n +i) \ (for fixed n — 5) in the Kondo lattice model, (a) Se- 
ries of results under usual OBC and deformed OBC by varying the 
system size L. (b) The variation of results of the deformed OBC 
(L = 44) for several choices of m, to be compared with that of the 
usual OBC with L — 40 and m = 600. In both panels, the result of 
PBC with L — 64 (m = 600 for non-Abelian Hilbert space) is given 
for comparison. 



to Eq. ([8j. Figure|6la) shows the comparison of the site depen- 
dence of electron density, (n.;), between usual OBC and de- 
formed OBC. In usual OBC, the electron density significantly 
deviates from 0.5 near the system boundary, and as a result, 
the density around the center slightly exceeds quarter filling. 
Such deviation due to the boundary effect is suppressed in the 
deformed case, and the electron density becomes nearly site 
independent. Thus, the results of the latter are guaranteed to 
be less affected by the missing bond between the two edges. 

The comparison of the local quantities of two boundary 
conditions is given in Fig. [7] The oscillation amplitude of 
electron hopping, (cjcj+i + h.c), and the nearest neighbor 
correlation of localized spins, (SfSf +1 ), are both suppressed 
by the deformation of interaction. It is interesting to find that 
the mean value of (S*S* +1 ) remains a positive finite value 
even under the deformation. In usual OBC, (SfSf +1 ) has 
a large twofold oscillation between the comparable negative 
and positive values, and the local spin correlations are char- 
acterized by classical q = it/ 2 antiferromagnetic state, ttll- 
However, by the comparison of the two results, one can con- 
clude that the antiferromagnetic correlation is only stabilized 
by the boundary effect (the missing bond in usual OBC). The 
remaining ferromagnetic local correlation seems to be more 



intrinsic than the antiferromagnetic one in the J/t — 1 KLM. 
In fact, the phase diagram given in Refl30l indicates that sys- 
tem at this parameter is very close to the ferromagnetic phase 
transition. 

Finally, we show in Figs. [8{a) and [HJb) the L and m 
dependencies of the two-point dimer correlation function, 
\(O n O„+i)\. The presented results show the following two 
important features: The decay of \ (O n O n +i) \ is always more 
rapid for larger system size L irrespective of whether we have 
usual OBC or deformed OBC. Having larger m also leads to 
the rapid decay in the deformed OBC. Let us discuss the im- 
plication of these two features. In general, the accuracy of 
DMRG is guaranteed by asymptotically restricting the num- 
ber of bases m needed to reproduce the quantum states. How- 
ever, this technical advantage may sometimes bring in a by- 
product when m is not taken large enough: an artificial sup- 
pression of intrinsic long-range quantum fluctuation which is 
required to suppress the boundary-induced symmetry braking. 
If we focus on the value of | (O n O n +i) | at the particular length 
i in Fig. [HI a), we always find smaller values for larger L. Since 
the calculation with larger L includes the component of quan- 
tum fluctuation with larger length scale, this fact indicates that 
the dimer correlation should indeed be suppressed if the sys- 
tem includes a proper long-range quantum fluctuation present 
in the bulk limit. The same discussion holds for the m de- 
pendence: by including a larger number of m, the more the 
fluctuation effect is included. If such fluctuation is intrinsic, 
the dimer correlation will be suppressed by the increase of m, 
which is in fact the case. 

Let us briefly comment on the relation of the above L and 
m scaling with the entanglement entropy. In the systems 
with small L and m, the number of eigenstates is limited to 
a small number. This smallness generally enhances the sep- 
aration of discrete energy levels and suppresses long-range 
quantum fluctuations, which severely influences the accuracy 
of long-range two-point correlation functions near the critical 
point. Meanwhile, the long-range quantum fluctuation gener- 
ates quantum entanglement between two regions in the sys- 
tem, so that the increase in L and m enhances the entangle- 
ment entropy. Therefore, not only the two-point correlation 
function but also the entanglement entropy shall serve as a 
measure of to what extent the long range quantum fluctuation 
is taken into account near the critical point. The detailed and 
quantitative analysis on the relation of entanglement entropy 
with the boundary effect is left as a future problem. 

We finally note that the results of these L and m scaling 
asymptotically approach the one under the PBC in the same 
panel of Fig. [8J The \(O n O n +i)\ of PBC clearly shows an 
algebraic decay characteristic of the Tomonaga-Luttinger liq- 
uid. As discussed in Ref. |Tl] (see Sec. II), the sin 2 deforma- 
tion is one of the ways to recover the translational invariance, 
which is demonstrated in Sec. Ill for the J1-J2 model. To re- 
cover perfectly the translational symmetry by the deformation 
in KLM at J = 1 is rather out of scope in the present calcula- 
tion, because the on-site Kondo singlet correlation competes 
with the inter-site hopping, and the applicability of the sin 2 
deformation is not clear. However, even in such a difficult 
case, the analysis of the asymptotic behavior of \(O n O n +i)\ 
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gives us the strong indication that the true bulk property of 
the system approaches the one given under PBC, and that the 
dimer correlation is only induced by the boundary effect. The 
good coincidence of \(O n O n +i) \ in Fig. [5] and Fig. [8] under 
different series of boundary conditions is regarded as its col- 
lateral evidence. 



V. SUMMARY AND DISCUSSION 

To summarize, we made a case study analysis on the den- 
sity matrix renormalization group calculation and proposed 
the systematic treatment to determine whether the symme- 
try breaking long range order exists or not. Since the two- 
point correlation function in the bulk limit gives the lower 
bound of the symmetry-breaking long-range order parameter 
(see Sec. Ill Al l, we focused on the two-point correlation func- 
tions and systematically studied their boundary condition de- 
pendence. We demonstrated the calculations on the two con- 
trasting cases, the J1-J2 model which has dimer long-range 
order, and the Kondo lattice model at quarter filling whose 
ground state has not reached a consensus yet due to numerical 
difficulty. 

Four different types of boundaries are adopted and their ef- 
fect on the expectation values of order operators are analyzed. 
In the usual OBC we always find oscillation of local quanti- 
ties, e.g., the nearest neighbor spin-spin correlation (SfSf +1 ) 
or bond kinetic energy (c|c,+i + h.c.) oscillates with the 
largest amplitude. The modified OBC relaxes the oscillation 
to some extent under the variation of imbalance between po- 
tentials on both edges. Then, the deformation of interactions 
and PBC basically recover the homogeneity of local quanti- 
ties. As the boundary effect fades out in such a way, the two- 
point correlation function also shows a systematic boundary 
condition dependence, and thus the degree of inhomogene- 
ity of the local quantities gives a measure of to what extent 
the boundary effect of OBC is included in the wave function. 
Therefore, by examining the correspondence of local quanti- 
ties and the two-point correlation function, one could analyze 
the overall features of the boundary effect as summarized in 



(III) of Sec. UFA] 

Such boundary condition dependence is particularly dis- 
tinct when the system is critical. When the system has 
symmetry-breaking long-range order (gapped), the two-point 
correlation function is insensitive to the boundary conditions. 
This fact is confirmed in the 3\- J2 model, where the two-point 
correlation functions reach the same value after the finite-size 
scaling even when the corresponding local correlations take 
the significantly different value with each other. By sharp 
contrast, when the system lacks a distinct gap, the boundary 
condition severely affects the two-point correlation function 
as we find in the Kondo lattice model:— In OBC it seems to 
behave convex upward and asymptotically approach a finite 
value with increasing distance, which is easily transformed to 
the algebraic decay, once either of the modified OBC, sin 2 de- 
formation, or PBC is adopted. The boundary condition depen- 
dence of the two-point correlation function (Figs(5]and[8]l is in 
good correspondence with that of the local quantities (Figs. [4] 
and [7]i, and such boundary-dependent behavior is clearly dif- 
ferent from that observed in the J\- J2 model. These results 
show that the dimer correlation observed in OBC is not an 
intrinsic property of the KLM and that it is difficult to con- 
clude that the symmetry-breaking long-range order does exist 
at quarter filling at J/t = 1. This fragile example which re- 
quires a difficult numerical treatment, provides us with a good 
lesson that numerical study, although it sometimes becomes 
as rigorous as analytical techniques, may have much analysis 
dependence, and one quite often needs care to understand the 
intrinsic physics lying behind. 

We conclude that the DMRG calculation, which is often 
performed on OBC, sometimes reaches misleading conclu- 
sions. While OBC is very useful to analyze the detailed nature 
of the well-known ground state, we argue that to have a reli- 
able conclusion on the determination of the ground state itself, 
the examination of the boundary condition dependence on the 
two-point correlation should be performed. 
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